Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1404 1616 (N39W39 L=339)   315     73        0 (1385/1583-1423/1647)  51.3 W36-W42 N36-N41
      2 1592 1580 (N34W55 L=355)  1247    195        0 (1504/1537-1681/1621) -21.8 W47-W64 N30-N37
      3 1060 1565 (N40W09 L=308)   325     46        1 (1041/1536-1077/1590)  53.1 W07-W10 N37-N41
      4 1716 1366 (N18W56 L=356)  1035    253        3 (1698/1233-1745/1462)  52.4 W50-W64 N10-N23
      5 1712  644 (S27W49 L=349)   356     31        2 (1689/ 626-1735/ 660)   0.0 W46-W52 S28-S25
      6 1267  574 (S26W12 L=312)   227     60        1 (1234/ 549-1288/ 594) -33.7 W10-W13 S27-S24
      7 1051  464 (S32E04 L=296)   632    231        1 (1007/ 368-1116/ 546)  53.1 E08-W01 S38-S26
      8 1254  473 (S33W11 L=311)   358     75        1 (1233/ 434-1270/ 518)  81.9 W09-W12 S36-S29
      9  937  284 (S46E17 L=283)   340     57        1 ( 902/ 261- 961/ 306) -59.0 E18-E14 S48-S43
  

    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