Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1302 1733 (N40W47 L=351)   339     56        0 (1270/1720-1337/1749) -45.0 W44-W49 N38-N41
      2 1047 1623 (N42W21 L=324)   259     58        1 (1019/1607-1078/1644) -33.7 W18-W22 N40-N43
      3  854 1453 (N36W01 L=304)   872    243        2 ( 759/1411- 974/1473) -24.8 E05-W08 N32-N39
      4 1772 1391 (N06W62 L=  5)   434    101        0 (1757/1335-1785/1431)  39.8 W58-W64 N03-N08
      5 1818 1127 (S09W55 L=359)   469     70        0 (1801/1095-1833/1168)  53.1 W52-W57 S10-S06
      6  229 1022 (N23E57 L=247)  1198    285        3 ( 190/ 919- 292/1146) -53.7 E62-E51 N15-N30
      7 1722  999 (S13W43 L=347)  1367    258        2 (1690/ 889-1746/1137)  61.7 W39-W47 S18-S05
      8 1857 1060 (S14W58 L=  2)   459     52        1 (1842/1025-1876/1088)  71.6 W56-W60 S15-S12
      9  798 1072 (N15E12 L=292)   187     94        0 ( 751/1051- 847/1087) -38.7 E14-E09 N12-N16
     10 1143  865 (S04W03 L=306)   571    210        0 (1118/ 775-1176/ 961)  56.3 W00-W06 S09-N00
     11  220  834 (N12E61 L=242)   868     92        5 ( 185/ 784- 258/ 870)  26.6 E66-E56 N08-N14
     12  283  806 (N09E55 L=249)   656     87        3 ( 254/ 764- 309/ 855)  39.8 E58-E51 N07-N12
     13  611  716 (S01E32 L=272)   247     72        0 ( 600/ 679- 628/ 753)  56.3 E33-E30 S03-N00
     14 1406  523 (S33W11 L=315)  1771    368        1 (1274/ 416-1508/ 654)  27.8 W00-W20 S36-S26
     15  302  457 (S13E73 L=230)  1431    129        2 ( 274/ 405- 326/ 509)  28.1 E82-E64 S17-S09
     16  858  227 (S39E39 L=265)  1688    272        1 ( 763/ 159- 977/ 304)  12.7 E56-E25 S43-S35
  

    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