Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1136 1648 (N43E04 L=319)   467     82        1 (1108/1608-1160/1686) -71.6 E06-E01 N39-N45
      2  361 1450 (N15E56 L=266)  1332    147        2 ( 333/1382- 391/1519) -66.0 E59-E53 N11-N20
      3  565 1480 (N21E40 L=283)   375     87        1 ( 529/1451- 606/1503)  21.8 E43-E37 N18-N22
      4 1684 1346 (N31W44 L=  7)   828    221        1 (1609/1267-1740/1427) -36.0 W37-W49 N26-N35
      5  870 1392 (N20E16 L=307)   294    100        0 ( 833/1364- 915/1424)  51.3 E17-E13 N17-N22
      6 1046 1274 (N15E03 L=320)   277    112        1 (1016/1234-1082/1325) -45.0 E05-W00 N13-N18
      7 1884 1250 (N28W71 L= 33)   711     53        1 (1862/1216-1901/1272)  29.7 W66-W76 N25-N29
      8 1800 1241 (N26W56 L= 19)   599     65        1 (1779/1206-1824/1269)  53.1 W53-W59 N23-N27
      9 1777 1120 (N17W53 L= 15)   398     61        2 (1747/1098-1812/1143) -18.4 W49-W56 N16-N18
     10 1204  856 (S07W14 L=336)   214    127        0 (1145/ 833-1257/ 876)   0.0 W09-W17 S08-S06
     11 1216  725 (S15W17 L=340)   844    265        0 (1145/ 641-1267/ 777)  90.0 W11-W19 S20-S11
     12 1431  664 (S16W33 L=355)   455    149        0 (1371/ 619-1477/ 708)  54.5 W29-W35 S19-S12
     13  593  636 (S31E24 L=299)  1226    367        2 ( 450/ 584- 758/ 699)  21.8 E36-E11 S37-S25
     14 1141  632 (S22W14 L=337)   213     35        1 (1124/ 610-1160/ 658) -45.0 W12-W15 S23-S20
     15 1554  444 (S28W53 L= 16)  3721    493        5 (1358/ 371-1740/ 560)  32.2 W38-W66 S34-S17
     16  994  379 (S43W11 L=334)  1417    273        4 ( 875/ 348-1125/ 418)  17.5 W01-W20 S47-S38
  

    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