Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  975 1705 (N32E26 L=181)   953    207        2 ( 913/1624-1031/1768) -23.2 E32-E18 N27-N36
      2  465 1680 (N18E62 L=146)  1758    143        7 ( 422/1615- 515/1739) -39.3 E71-E54 N13-N22
      3  683 1603 (N17E39 L=169)   753    227        2 ( 596/1568- 777/1635)  22.6 E45-E32 N15-N21
      4  584 1551 (N12E43 L=164)   818    196        0 ( 537/1510- 664/1629) -14.9 E50-E35 N10-N16
      5 1404 1388 (N25W12 L=219)   563    224        0 (1373/1295-1422/1488) -77.0 W08-W13 N18-N31
      6 1603 1367 (N31W27 L=234)   464    101        0 (1575/1316-1623/1411) -45.0 W23-W28 N27-N32
      7  899 1375 (N09E17 L=190)   299    111        0 ( 851/1342- 947/1401)  -8.1 E20-E13 N08-N10
      8  794 1352 (N05E22 L=185)   207     66        0 ( 771/1325- 813/1385)  90.0 E22-E21 N03-N07
      9  278 1322 (S08E59 L=148)   710    135        0 ( 246/1257- 310/1384)  82.9 E60-E57 S11-S03
     10  358 1286 (S08E50 L=158)   333     88        0 ( 333/1248- 380/1326)  78.7 E50-E48 S10-S05
     11 1195  860 (S11W14 L=221)   187     42        0 (1166/ 842-1214/ 881)  56.3 W12-W15 S12-S09
     12  825  758 (S27E04 L=203)   637    160        1 ( 772/ 725- 897/ 798)   0.0 E08-W00 S29-S25
     13 1379  420 (S29W45 L=252)   618    125        1 (1328/ 402-1449/ 451)  16.7 W40-W50 S30-S26
  

    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