Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  664 1770 (N33E48 L=  0)   562     73        1 ( 632/1746- 703/1795) -11.3 E52-E42 N31-N34
      2 1241 1756 (N46E04 L= 44)   532     83        0 (1209/1728-1277/1784)  51.3 E05-W00 N42-N48
      3  439 1665 (N21E60 L=348)  2474    188        8 ( 342/1618- 503/1736)  38.7 E72-E52 N18-N27
      4  922 1638 (N28E21 L= 26)   460    105        2 ( 890/1582- 953/1682) -45.0 E23-E17 N24-N31
      5 1502 1468 (N31W21 L= 69)  1894    477        2 (1318/1311-1643/1596) -20.3 W06-W33 N24-N35
      6 1120 1400 (N17E02 L= 45)   331    122        0 (1086/1353-1163/1445) -33.7 E05-W01 N14-N18
      7  661 1244 (S01E26 L= 22)   297     65        0 ( 639/1216- 680/1278)  71.6 E27-E24 S03-N00
      8 1811 1085 (N16W51 L= 99)  3237    439        5 (1731/ 882-1892/1257) -55.0 W41-W61 N04-N24
      9  836 1163 (S03E14 L= 34)   202     47        2 ( 808/1145- 865/1184)  14.0 E15-E11 S04-S01
     10  496  929 (S22E32 L= 16)   379     82        1 ( 484/ 883- 510/ 977) -71.6 E33-E30 S25-S19
     11  161  838 (S32E68 L=340)  1397    195        2 ( 129/ 742- 211/ 914)  26.6 E81-E58 S38-S27
     12 1709  648 (S11W55 L=103)   380     47        0 (1691/ 629-1732/ 675) -26.6 W53-W58 S12-S09
     13  517  629 (S41E28 L= 20)   303     72        0 ( 483/ 608- 559/ 644)  29.7 E31-E24 S43-S39
     14 1477  493 (S26W45 L= 93)  1739    332        4 (1325/ 414-1590/ 580)   0.0 W30-W57 S28-S23
     15 1457  342 (S35W55 L=103)  1066    128        1 (1409/ 309-1512/ 408)  -9.5 W45-W63 S37-S32
     16 1376  281 (S41W55 L=103)   548     62        1 (1353/ 250-1398/ 314) -18.4 W50-W61 S42-S39
  

    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