Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  866 1871 (N42E48 L=148)  1052    134        2 ( 816/1829- 928/1910)  -8.5 E59-E38 N40-N45
      2 1114 1644 (N32E14 L=182)  1174    216        2 (1050/1552-1168/1730) -31.0 E22-E07 N27-N36
      3  553 1685 (N20E54 L=143)   557    115        1 ( 525/1643- 585/1721) -63.4 E56-E49 N16-N22
      4  496 1686 (N19E59 L=137)   428     48        2 ( 475/1667- 523/1709)   0.0 E62-E56 N17-N20
      5  770 1566 (N17E31 L=165)   521    169        1 ( 727/1525- 841/1619) -21.8 E35-E25 N15-N20
      6  702 1504 (N12E33 L=163)   479    139        1 ( 660/1457- 754/1557) -21.8 E37-E27 N10-N14
      7 1240 1381 (N19W02 L=199)   336     98        1 (1215/1343-1276/1431) -45.0 W00-W05 N17-N21
      8  397 1335 (S05E48 L=148)   858    205        1 ( 367/1239- 424/1425) -68.2 E51-E45 S10-N00
      9 1544 1308 (N25W23 L=220)   795    223        0 (1509/1223-1565/1411) -61.4 W18-W25 N18-N30
     10 1699 1317 (N31W36 L=232)   875    205        0 (1626/1247-1746/1389) -31.0 W28-W40 N27-N33
     11  498 1269 (S06E38 L=158)   716    217        0 ( 442/1183- 549/1361)  80.5 E39-E36 S12-N00
     12  967 1281 (N06E10 L=186)   436    125        0 ( 930/1226-1005/1323)  82.9 E11-E09 N01-N09
     13 1571 1202 (N19W27 L=224)   240     29        0 (1552/1180-1590/1222)  71.6 W25-W28 N17-N20
     14 1552  922 (N03W33 L=230)   222     43        2 (1528/ 904-1579/ 944) -14.0 W31-W35 N01-N04
     15  182  870 (S36E64 L=132)   986    164        1 ( 133/ 829- 237/ 909)  17.5 E76-E55 S39-S33
     16 1345  778 (S11W25 L=222)   209     20        0 (1322/ 761-1364/ 797)  18.4 W23-W26 S11-S09
     17  968  685 (S27W07 L=203)   448    150        1 ( 898/ 661-1029/ 720)  12.5 W01-W10 S28-S25
     18 1463  358 (S29W56 L=253)   943    123        1 (1412/ 336-1531/ 391)  20.0 W50-W63 S30-S26
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y