Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1229 1640 (N34W33 L=136)   584    142        1 (1159/1611-1286/1661) -11.3 W27-W37 N30-N35
      2  864 1395 (N29W00 L=104)   527    117        0 ( 804/1373- 928/1416) -33.7 E03-W03 N26-N31
      3 1585 1335 (N07W42 L=146)   294     82        0 (1565/1291-1607/1366)  31.0 W39-W44 N04-N08
      4 1726 1265 (S01W51 L=154)  1034    252        1 (1642/1156-1800/1363) -66.8 W47-W54 S08-N06
      5  859 1259 (N21E04 L=100)   261     71        1 ( 837/1220- 893/1296)  31.0 E05-W00 N19-N22
      6 1031 1189 (N13W04 L=108)   178     56        0 ( 998/1170-1065/1207) -45.0 W02-W06 N11-N14
      7 1556 1067 (S07W32 L=136)   896    240        2 (1519/ 946-1580/1182)  72.9 W29-W34 S13-N00
      8  200 1023 (N22E57 L= 47)  2090    479        3 (  92/ 985- 310/1089) -20.0 E68-E46 N18-N27
      9  796 1066 (N12E12 L= 92)   236     86        0 ( 750/1049- 842/1080)  -9.5 E14-E08 N11-N12
     10  403  958 (N15E40 L= 64)   434    103        1 ( 360/ 905- 432/ 994) -63.4 E42-E38 N11-N17
     11 1794  946 (S21W48 L=152)   341     84        0 (1764/ 923-1824/ 963)   0.0 W44-W51 S21-S19
     12 1216  801 (S13W05 L=109)   323     89        1 (1165/ 774-1262/ 825)   0.0 W01-W08 S13-S11
     13 1394  584 (S31W11 L=115)   541    114        1 (1365/ 539-1420/ 643)  74.1 W09-W12 S34-S27
     14  358  570 (S07E57 L= 47)  1078    196        1 ( 289/ 535- 442/ 612) -26.6 E62-E49 S09-S03
     15  940  498 (S24E19 L= 84)   311    124        0 ( 912/ 445- 962/ 552)  39.8 E22-E16 S26-S20
     16 1115  467 (S31E10 L= 94)   245     69        0 (1090/ 431-1138/ 498)  31.0 E12-E07 S32-S29
  

    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