Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1363 1629 (N42W33 L=289)   389     83        0 (1321/1611-1400/1643) -26.6 W29-W36 N40-N43
      2 1230 1605 (N41W20 L=277)   275     58        0 (1210/1576-1248/1629)  45.0 W18-W22 N39-N43
      3 1759 1279 (N15W56 L=313)   741     97        0 (1718/1251-1801/1304) -23.2 W52-W60 N13-N16
      4  309 1164 (N14E51 L=205)   328     58        0 ( 277/1148- 335/1180)  -9.5 E54-E48 N13-N15
      5 1062  869 (S05W02 L=258)   283    115        0 (1003/ 844-1114/ 887) -21.8 E01-W04 S07-S03
      6  405  846 (S05E43 L=214)   283     59        0 ( 390/ 816- 424/ 879) -90.0 E44-E41 S07-S03
      7 1420  642 (S22W25 L=281)   485     53        1 (1398/ 613-1450/ 672)  45.0 W23-W27 S23-S19
      8  557  600 (S21E35 L=222)   331     59        0 ( 531/ 572- 580/ 627) -45.0 E36-E33 S23-S19
      9 1292  529 (S29W17 L=273)   737    142        2 (1233/ 482-1350/ 577)  36.9 W12-W20 S32-S25
     10  539  515 (S27E39 L=217)   333    110        0 ( 516/ 472- 577/ 554)  40.6 E41-E34 S30-S24
     11  931  498 (S29E09 L=247)   713    255        0 ( 839/ 444-1045/ 552)  20.6 E16-W00 S33-S25
     12 1566  392 (S43W47 L=303)   663     49        0 (1540/ 369-1593/ 420) -26.6 W42-W52 S45-S39
  

    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