Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  453 1514 (N36E45 L=349)  1659    558        2 ( 325/1436- 581/1607) -18.4 E62-E32 N30-N42
      2 1656 1533 (N34W53 L= 87)   714     79        1 (1617/1502-1685/1563) -45.0 W50-W55 N31-N36
      3  334 1288 (N20E50 L=344)  2894    685        6 ( 214/1089- 422/1442) -37.9 E67-E40 N08-N29
      4  876 1386 (N27E09 L= 25)   317     78        1 ( 833/1371- 921/1403)   9.5 E12-E06 N26-N28
      5 1836 1300 (N17W65 L= 99)   647    114        0 (1803/1259-1862/1338) -45.0 W61-W67 N14-N19
      6 1079 1194 (N14W04 L= 38)   351    138        1 (1027/1132-1142/1249) -45.0 W00-W07 N10-N17
      7 1487  982 (S00W30 L= 64)   272     93        1 (1442/ 955-1530/1007)   9.5 W26-W32 S01-N01
      8  440  914 (S03E39 L=355)   280     81        1 ( 417/ 888- 470/ 944) -63.4 E40-E36 S04-N00
      9  726  828 (S08E19 L= 15)   553    154        0 ( 652/ 783- 796/ 863) -24.0 E23-E14 S10-S05
     10  338  693 (S17E51 L=344)  1035    250        1 ( 301/ 607- 377/ 796)  50.2 E56-E45 S22-S10
     11 1373  578 (S26W23 L= 58)   609    143        2 (1306/ 560-1426/ 608)  14.0 W18-W27 S26-S23
     12  550  539 (S27E35 L=359)   874    233        1 ( 487/ 503- 646/ 597)  25.0 E41-E26 S29-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