Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1546 1500 (N40W29 L=148)   327     63        0 (1525/1476-1574/1517)  36.9 W27-W31 N37-N41
      2 1239 1390 (N25W07 L=125)   625    145        1 (1183/1331-1312/1446) -26.6 W02-W12 N22-N27
      3 1765 1399 (N37W50 L=169)  1033    193        0 (1711/1356-1821/1432) -16.7 W44-W56 N35-N39
      4 1856 1311 (N33W59 L=178)   709    178        1 (1823/1294-1891/1324)  12.5 W55-W64 N31-N34
      5 1700 1079 (N15W44 L=163)  1491    297        4 (1650/ 967-1781/1161)  77.0 W40-W50 N07-N20
      6 1485  842 (S03W32 L=151)   523    165        1 (1407/ 818-1557/ 873)  26.6 W26-W36 S05-N00
      7  195  777 (S31E59 L= 60)   607     76        0 ( 168/ 745- 225/ 805) -36.9 E62-E55 S33-S29
      8 1802  601 (S11W70 L=188)   757    243        0 (1769/ 581-1838/ 629)  -4.1 W63-W77 S12-S10
      9 1474  568 (S20W39 L=158)   332     77        0 (1459/ 538-1489/ 603) -53.1 W37-W41 S21-S17
  

    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