Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  821 1799 (N40E36 L=328)   828    107        2 ( 768/1771- 872/1822)  -5.2 E41-E30 N37-N41
      2 1162 1703 (N39E06 L=358)  1016    143        3 (1113/1628-1194/1756) -52.1 E10-E02 N34-N43
      3 1061 1725 (N39E14 L=350)   480    114        0 (1008/1692-1102/1753) -12.5 E18-E09 N36-N40
      4 1020 1360 (N13E07 L=357)   594    290        1 ( 899/1314-1140/1395)  17.1 E13-W00 N08-N16
      5  207 1155 (S11E61 L=303)   439     56        2 ( 189/1124- 221/1188) -53.1 E63-E58 S12-S08
      6  299  981 (S20E49 L=314)   320     83        0 ( 289/ 934- 308/1030) -90.0 E50-E48 S23-S17
      7  261  776 (S33E55 L=309)   721    124        0 ( 219/ 732- 296/ 835)  32.5 E62-E50 S36-S29
      8 1287  639 (S24W25 L= 29)   787    256        0 (1260/ 528-1326/ 739) -69.0 W20-W28 S30-S17
      9 1595  678 (S14W45 L= 48)   296     74        2 (1579/ 638-1611/ 719) -68.2 W42-W46 S16-S11
     10 1693  666 (S12W53 L= 57)   534     84        3 (1665/ 629-1728/ 702) -23.2 W50-W57 S13-S10
     11 1674  596 (S17W55 L= 59)   527      0       -1 (1651/ 540-1704/ 640) -66.8 W52-W57 S20-S13
     12 1474  442 (S31W48 L= 51)   918    170        1 (1411/ 398-1548/ 503)  66.0 W43-W51 S34-S25
     13 1587  396 (S30W62 L= 66)  1030    166        1 (1516/ 357-1653/ 443)  56.3 W56-W65 S34-S25
  

    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