Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1087 1685 (N48W21 L=327)   489     92        0 (1058/1649-1133/1706)   7.1 W16-W24 N45-N49
      2 1196 1591 (N39W25 L=331)   777    121        2 (1160/1545-1244/1641) -77.5 W21-W26 N34-N43
      3 1583 1523 (N26W52 L=359)   455     98        0 (1540/1493-1624/1554) -51.3 W49-W54 N23-N28
      4 1429 1381 (N21W33 L=340)   225     35        0 (1401/1366-1453/1393) -33.7 W31-W35 N19-N21
      5  173 1149 (N23E68 L=238)   898    152        1 ( 149/1067- 189/1210) -52.1 E74-E64 N17-N26
      6  558 1086 (N16E29 L=277)  1023    391        1 ( 415/1002- 713/1168) -29.7 E40-E19 N09-N22
      7 1731  865 (S16W46 L=353)   594    205        0 (1651/ 798-1789/ 938)  24.8 W38-W52 S18-S12
      8 1516  864 (S11W29 L=335)   796    150        2 (1476/ 798-1565/ 932)  40.6 W25-W32 S14-S07
      9 1811  841 (S20W55 L=  1)  1357    172        2 (1791/ 756-1827/ 928) -76.0 W52-W57 S25-S13
     10 1012  795 (S07E04 L=302)   802    196        2 ( 974/ 698-1058/ 877)  59.0 E07-W00 S12-S02
     11 1508  647 (S25W27 L=333)   257     97        0 (1481/ 611-1531/ 694)  51.3 W24-W28 S27-S22
  

    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