Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1059 1717 (N36E15 L= 88)   711    126        2 ( 994/1679-1115/1750) -10.3 E21-E10 N34-N38
      2  423 1561 (N14E53 L= 51)   329     78        1 ( 391/1543- 455/1579)  26.6 E56-E50 N12-N15
      3 1340 1509 (N29W09 L=113)   740    265        2 (1232/1459-1456/1543)  28.1 W02-W17 N22-N31
      4 1589 1500 (N35W28 L=132)   958    212        0 (1507/1437-1665/1533)   4.1 W20-W34 N32-N37
      5 1121 1248 (N07W01 L=105)   283     83        0 (1079/1228-1167/1265)  36.9 E01-W03 N05-N08
      6  215 1217 (S09E60 L= 44)   503    117        0 ( 187/1161- 240/1263)  71.6 E62-E58 S12-S06
      7  344 1133 (S13E46 L= 58)   467    110        1 ( 292/1104- 401/1161)  -5.7 E50-E40 S13-S11
      8  618 1096 (S11E25 L= 78)   236    104        0 ( 594/1044- 635/1143)  90.0 E26-E24 S14-S07
      9 1848 1086 (N17W54 L=158)   897    107        0 (1803/1043-1888/1130) -11.3 W49-W59 N14-N18
     10  926  961 (S13E04 L= 99)   849    497        0 ( 718/ 871-1115/1039)  36.9 E15-W05 S21-S05
     11 1873  876 (N06W62 L=166)  1054    216        1 (1849/ 768-1890/ 969) -63.4 W58-W66 N00-N10
     12  324  946 (S24E46 L= 58)   313     58        0 ( 289/ 931- 356/ 963)  18.4 E49-E42 S25-S22
     13 1426  863 (S07W27 L=131)   206     56        0 (1393/ 846-1453/ 878)   0.0 W24-W28 S07-S06
     14 1162  824 (S15W12 L=116)   204     87        1 (1140/ 782-1184/ 856)  90.0 W11-W13 S18-S13
     15 1191  669 (S23W18 L=122)   297    126        0 (1132/ 655-1242/ 689)  15.9 W14-W21 S24-S21
     16 1805  592 (S11W69 L=172)   617     96        4 (1775/ 567-1828/ 623) -15.9 W64-W72 S12-S09
     17 1180  584 (S29W20 L=124)   271     77        0 (1143/ 568-1219/ 603)  45.0 W18-W22 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