Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1445 1750 (N37W61 L=289)  1549    225        1 (1351/1711-1543/1779) -12.3 W48-W72 N33-N40
      2 1566 1363 (N13W43 L=271)  1242    294        4 (1525/1224-1631/1490)  66.4 W39-W48 N03-N19
      3 1062 1088 (N10W04 L=232)   208    100        0 (1034/1039-1081/1131) -81.9 W02-W04 N06-N13
      4 1374 1087 (N03W22 L=250)   214     42        2 (1349/1061-1392/1113)  71.6 W20-W22 N01-N04
      5 1213  976 (N00W10 L=238)   931    314        0 (1143/ 850-1325/1035)  19.7 W03-W17 S05-N03
      6  225  923 (N16E59 L=169)   644    109        1 ( 187/ 880- 260/ 957) -45.0 E63-E55 N12-N17
      7 1571  891 (S14W30 L=258)   430    126        0 (1541/ 838-1621/ 926)  26.6 W27-W34 S16-S11
      8 1586  814 (S19W30 L=258)   400     65        0 (1559/ 779-1607/ 845)  36.9 W27-W31 S20-S16
      9  972  793 (S05E08 L=220)   187     63        0 ( 948/ 769-1000/ 820) -63.4 E08-E06 S06-S02
     10  241  723 (N04E63 L=166)   593    106        0 ( 186/ 712- 292/ 736)  -4.1 E70-E56 N02-N04
     11  689  651 (S07E28 L=200)   355     89        1 ( 674/ 610- 709/ 701)  68.2 E29-E26 S09-S04
     12  407  491 (S13E56 L=172)  8321    555       17 ( 318/ 340- 506/ 666)  50.4 E78-E46 S23-N00
     13 1830  628 (S41W58 L=286)   790     85        1 (1800/ 598-1864/ 655)  31.0 W51-W64 S43-S37
     14  498  533 (S11E46 L=182)   309     23        1 ( 481/ 511- 518/ 556)  56.3 E48-E44 S12-S09
     15 1264  476 (S32W02 L=230)   251     86        0 (1247/ 446-1294/ 513) -78.7 W00-W03 S34-S29
     16  636  358 (S26E45 L=183)  1719    353        3 ( 537/ 244- 761/ 495)  21.5 E67-E29 S33-S18
  

    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