Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1304 1696 (N35W44 L=148)  1406    255        1 (1207/1661-1424/1720) -29.7 W36-W52 N29-N38
      2  390 1409 (N45E37 L= 67)  1553    201        4 ( 303/1344- 480/1461) -21.4 E50-E27 N39-N51
      3 1852 1247 (S07W63 L=167)   476     54        2 (1835/1213-1870/1270)  71.6 W60-W66 S08-S05
      4 1553 1159 (S02W34 L=138)   538    132        0 (1529/1100-1577/1219) -90.0 W32-W35 S05-N02
      5 1820 1039 (S17W53 L=157)  1247    284        0 (1740/ 941-1892/1133)  12.8 W42-W64 S19-S14
      6  156 1066 (N28E66 L= 38)  1715    150        3 (  92/1006- 217/1106) -17.7 E81-E57 N23-N31
      7 1604  676 (S31W27 L=131)  1369    307        2 (1560/ 555-1635/ 823)  73.6 W22-W29 S38-S21
      8  250  734 (N07E61 L= 43)   587     77        1 ( 227/ 708- 279/ 763)  15.9 E64-E56 N05-N09
      9  926  591 (S16E18 L= 87)   285     98        0 ( 916/ 542- 942/ 635)  56.3 E19-E15 S19-S13
     10  788  539 (S16E28 L= 77)   241     76        1 ( 753/ 507- 806/ 560) -76.0 E29-E26 S17-S13
     11  760  435 (S21E34 L= 70)   963    172        2 ( 698/ 376- 823/ 497)  17.1 E40-E27 S24-S18
     12 1359  411 (S41W03 L=107)   417    108        0 (1312/ 390-1412/ 428)  -7.1 E01-W07 S42-S40
  

    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