Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  672 1751 (N44E47 L=296)   531     70        3 ( 644/1716- 706/1771)  90.0 E51-E43 N41-N46
      2  453 1619 (N31E58 L=285)  1304    153        2 ( 404/1562- 490/1684) -71.6 E63-E51 N26-N35
      3 1242 1573 (N39W08 L=351)   262     67        0 (1211/1543-1268/1593) -36.9 W05-W10 N37-N40
      4 1121 1551 (N36E01 L=342)   315     95        0 (1087/1514-1154/1586) -38.7 E03-W02 N33-N38
      5 1432 1417 (N30W24 L=  6)   377    119        0 (1395/1359-1463/1457) -50.2 W20-W26 N26-N32
      6  525 1297 (N11E36 L=306)   571    162        1 ( 512/1222- 545/1363) -71.6 E38-E33 N06-N15
      7  918 1233 (N12E09 L=333)   424    207        0 ( 884/1137- 936/1313) -63.4 E12-E06 N06-N16
      8 1228 1172 (N12W11 L=353)   657    237        2 (1178/1076-1303/1262)  65.2 W08-W14 N05-N18
      9 1830 1066 (N13W59 L= 42)   698     80        1 (1789/1035-1863/1084)   6.3 W54-W63 N10-N13
     10  165 1033 (S10E67 L=276)  1003    116        3 ( 142/ 979- 194/1080)  45.0 E70-E62 S13-S06
     11 1134  896 (S06W08 L=351)   188     70        0 (1117/ 857-1147/ 927)  90.0 W06-W09 S07-S03
     12 1392  584 (S22W31 L= 13)  2539    577        7 (1168/ 480-1626/ 701)  31.8 W17-W46 S31-S12
     13  833  545 (S32E07 L=336)   475     97        1 ( 785/ 528- 893/ 569)   0.0 E10-E02 S33-S30
  

    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