Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  955 1773 (N43E19 L=166)   286     80        0 ( 932/1750- 987/1797)   0.0 E21-E16 N40-N44
      2  439 1530 (N19E49 L=137)  1113    187        4 ( 386/1462- 477/1600)  78.7 E52-E44 N13-N23
      3 1295 1464 (N25W11 L=196)   586    101        3 (1258/1430-1329/1510)  71.6 W08-W13 N22-N28
      4 1875 1265 (N23W62 L=248)  2183    293        2 (1757/1169-1938/1351) -21.8 W48-W73 N17-N26
      5 1360 1307 (N16W17 L=202)   288     94        1 (1318/1281-1401/1331)   0.0 W14-W19 N14-N16
      6  676 1228 (N02E24 L=162)   326    145        0 ( 619/1192- 731/1253)  -7.1 E27-E19 N00-N03
      7 1149 1161 (N04W06 L=191)   279    128        2 (1108/1100-1182/1202) -39.8 W02-W08 N00-N06
      8  796 1042 (S08E14 L=172)   224    114        1 ( 772/ 988- 820/1083)  90.0 E15-E12 S11-S05
      9  920  901 (S15E05 L=181)   196    150        0 ( 867/ 870- 948/ 958)  56.3 E07-E03 S16-S10
     10 1578  862 (S07W37 L=223)   436     89        1 (1541/ 835-1621/ 885)  18.4 W34-W40 S08-S05
     11  861  774 (S23E07 L=179)   568    200        2 ( 795/ 751- 951/ 813)   5.7 E11-E01 S24-S20
     12 1474  702 (S18W34 L=219)   633    267        2 (1348/ 653-1584/ 738)  29.7 W26-W41 S22-S14
  

    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