Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1536 1680 (N29W61 L=168)  1278    161        2 (1479/1623-1607/1716)   0.0 W51-W72 N26-N29
      2 1655 1509 (N16W58 L=164)   535    112        0 (1610/1480-1697/1535) -12.5 W52-W61 N14-N17
      3 1255 1448 (N25W26 L=133)   632    150        1 (1193/1397-1321/1489)  -6.3 W21-W30 N22-N26
      4  241 1274 (N39E57 L= 50)   695    152        2 ( 229/1204- 257/1360) -68.2 E59-E53 N34-N44
      5  992 1287 (N23W05 L=112)   529    172        1 ( 960/1209-1029/1352)  54.5 W02-W08 N17-N25
      6 1268 1257 (N14W20 L=127)   195    139        0 (1231/1204-1300/1312) -81.9 W19-W20 N10-N17
      7 1394 1204 (N08W26 L=133)   248    106        0 (1356/1164-1428/1249)  23.2 W22-W29 N06-N09
      8 1877 1126 (S12W62 L=169)  3698    261        9 (1810/1039-1934/1218)  33.7 W53-W75 S17-S06
      9 1693  777 (S25W37 L=144)  3709    848        7 (1545/ 508-1829/1088)  44.0 W21-W53 S40-S09
     10 1052  854 (S03E02 L=104)   190     77        0 (1015/ 839-1094/ 873) -36.9 E04-W00 S04-S01
     11 1758  752 (S30W43 L=150)   486    114        0 (1736/ 695-1777/ 800)  50.2 W40-W45 S32-S26
     12  915  664 (S11E15 L= 92)   236    118        0 ( 864/ 647- 973/ 690) -38.7 E17-E12 S12-S08
     13 1328  613 (S24W08 L=115)   220     89        0 (1294/ 580-1364/ 637)  18.4 W04-W10 S25-S23
     14 1187  492 (S28E04 L=103)   247     63        0 (1170/ 469-1215/ 526) -90.0 E05-E02 S30-S25
     15  365  443 (S14E65 L= 42)  2158    126        4 ( 318/ 383- 401/ 498)  48.4 E72-E58 S18-S09
     16  582  419 (S19E46 L= 61)   360     78        1 ( 568/ 375- 592/ 463)  50.2 E48-E43 S21-S15
  

    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