Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  451 1554 (N12E52 L=161)  1407    308        3 ( 347/1460- 540/1639)  63.4 E56-E48 N04-N18
      2  777 1577 (N19E29 L=183)   498    172        0 ( 718/1520- 828/1635) -24.4 E35-E24 N16-N22
      3  672 1522 (N14E34 L=178)   681    199        1 ( 638/1439- 710/1592) -73.3 E39-E31 N08-N18
      4 1143 1485 (N22E05 L=208)   283     68        0 (1116/1445-1166/1513) -45.0 E06-E02 N19-N23
      5 1763 1444 (N40W45 L=257)   747     92        1 (1728/1409-1797/1487) -33.7 W41-W49 N37-N42
      6  737 1363 (N05E25 L=188)   240     74        0 ( 700/1334- 774/1385)  -9.5 E27-E21 N04-N06
      7  493 1155 (S11E35 L=177)   461    172        0 ( 428/1101- 551/1190) -10.3 E40-E29 S13-S09
      8  170 1071 (S21E65 L=148)   631     78        2 ( 148/1028- 196/1107) -78.7 E67-E61 S23-S18
      9  634  994 (S18E22 L=190)   255     87        1 ( 611/ 960- 666/1035)  68.2 E23-E21 S19-S14
  

    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