Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  313 1602 (N17E67 L= 79)   779    165        1 ( 295/1571- 336/1641) -78.7 E69-E64 N15-N20
      2  367 1605 (N18E60 L= 86)   536     71        2 ( 337/1583- 393/1633)  45.0 E63-E57 N16-N20
      3 1093 1332 (N12E02 L=144)   809    254        1 (1069/1213-1123/1455) -66.8 E05-W02 N05-N19
      4 1718 1195 (N18W41 L=187)   436    124        1 (1705/1128-1727/1259) -76.0 W39-W42 N14-N22
      5 1358 1159 (N07W17 L=163)   381     96        0 (1325/1110-1390/1210)  74.1 W15-W18 N04-N11
      6  172 1163 (S11E64 L= 82)   458     74        0 ( 158/1120- 189/1197)  76.0 E66-E62 S13-S09
      7 1671  615 (S15W53 L=199)  1605    372        4 (1521/ 508-1790/ 771)  60.3 W45-W58 S24-S03
      8 1209  661 (S24W19 L=165)   862    226        0 (1083/ 638-1315/ 700)  28.3 W11-W24 S27-S19
      9 1421  537 (S26W37 L=184)   270     56        1 (1392/ 523-1450/ 555)  14.0 W35-W39 S27-S24
     10 1471  271 (S40W57 L=203)  2243    390        4 (1365/ 248-1564/ 299)  18.4 W47-W65 S43-S36
  

    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