Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  858 1801 (N37E39 L=102)  2354    399        3 ( 785/1705- 912/1890) -31.8 E53-E32 N29-N42
      2  431 1690 (N18E65 L= 77)   727    105        2 ( 390/1665- 465/1718)  45.0 E67-E62 N16-N20
      3  447 1565 (N11E54 L= 87)  1799    302        4 ( 352/1486- 539/1645)  63.4 E60-E48 N03-N16
      4 1160 1418 (N20E03 L=138)   839    162        2 (1115/1362-1222/1464) -20.6 E07-W01 N16-N22
      5 1437  623 (S15W37 L=178)  1113    224        0 (1314/ 578-1526/ 678)  50.2 W30-W40 S21-S09
      6 1440  559 (S19W40 L=182)   413     72        2 (1407/ 530-1481/ 585)  53.1 W38-W42 S20-S16
      7 1732  510 (S11W67 L=209)   705    116        4 (1705/ 487-1761/ 541)  63.4 W64-W71 S12-S08
      8 1541  430 (S22W54 L=196)  1204    297        2 (1469/ 373-1596/ 498) -10.6 W45-W61 S23-S20
      9 1611  465 (S18W58 L=200)   472     64        2 (1590/ 441-1634/ 498) -76.0 W55-W60 S19-S15
     10 1049  358 (S43W27 L=169)  1313    362        1 ( 894/ 329-1222/ 398)  23.5 W15-W39 S48-S38
  

    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