Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  643 1700 (N60E24 L= 24)   545     57        1 ( 616/1669- 668/1721)  68.2 E29-E20 N56-N61
      2  513 1651 (N57E41 L=  7)  2062    223        3 ( 428/1591- 603/1702)   4.6 E54-E28 N52-N61
      3  835 1560 (N44E04 L= 44)   581    117        2 ( 769/1535- 894/1578) -24.0 E09-W00 N41-N46
      4 1278 1533 (N33W29 L= 77)   255     31        1 (1253/1511-1299/1549) -33.7 W27-W30 N31-N34
      5 1039 1523 (N37W11 L= 59)   431    145        0 ( 973/1508-1108/1536)  -5.2 W05-W16 N35-N38
      6 1136 1302 (N21W12 L= 60)   662    181        2 (1087/1249-1212/1334)  12.5 W07-W17 N18-N23
      7  238 1270 (N31E61 L=347)   417     68        0 ( 207/1252- 277/1286)   5.7 E65-E55 N30-N32
      8  731 1065 (N13E17 L= 31)   873    266        2 ( 677/ 976- 759/1172)  76.0 E22-E14 N08-N20
      9  261 1087 (N20E55 L=353)  1032    137        3 ( 233/1023- 301/1134)  66.8 E58-E50 N15-N22
     10 1870  987 (S13W63 L=111)   511     70        0 (1850/ 948-1886/1016) -68.2 W60-W65 S15-S10
     11 1793  842 (S20W52 L=100)  1284    302        0 (1721/ 729-1851/ 967)  40.6 W44-W59 S25-S13
     12 1559  873 (S12W32 L= 80)   391    134        1 (1527/ 826-1592/ 943)  50.2 W28-W34 S14-S08
     13 1547  695 (S23W29 L= 77)  1031    205        2 (1459/ 645-1615/ 765)  17.1 W22-W35 S25-S19
     14 1098  672 (S16E01 L= 47)   604    149        2 (1068/ 603-1148/ 728)  50.2 E03-W02 S20-S12
     15 1537  577 (S31W28 L= 77)   311     68        0 (1507/ 554-1563/ 605)  21.8 W25-W30 S32-S29
     16 1183  521 (S28W02 L= 50)   630    123        1 (1154/ 476-1226/ 567) -66.8 W00-W03 S31-S24
  

    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