Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  755 1705 (N39E26 L=298)  1048    184        4 ( 699/1612- 802/1764) -45.0 E33-E20 N31-N43
      2  616 1711 (N39E39 L=286)   384     32        1 ( 595/1687- 637/1731)  90.0 E41-E36 N36-N40
      3  465 1576 (N28E46 L=279)   577     93        1 ( 419/1547- 513/1594)  -5.7 E50-E40 N26-N29
      4  318 1152 (N01E48 L=276)   508     77        0 ( 277/1120- 355/1186)  18.4 E51-E44 S01-N02
      5  500  845 (S17E33 L=291)   333     99        0 ( 486/ 793- 520/ 896)  71.6 E34-E31 S20-S14
      6  654  804 (S19E23 L=302)   251     96        0 ( 614/ 763- 697/ 836)  33.7 E25-E19 S22-S17
      7  314  749 (S23E51 L=274)   972    149        1 ( 266/ 687- 369/ 794)  28.3 E57-E44 S27-S20
      8 1710  635 (S22W54 L= 18)   317     23        1 (1686/ 620-1733/ 652)   0.0 W51-W56 S23-S21
      9  667  503 (S39E24 L=300)   874    124        2 ( 602/ 467- 725/ 545)  20.0 E30-E19 S41-S36
     10 1753  498 (S30W69 L= 33)  1046     68        3 (1718/ 466-1780/ 538)  68.2 W64-W76 S32-S27
     11 1265  414 (S42W24 L=349)   591    153        0 (1229/ 388-1328/ 468) -21.8 W19-W29 S43-S37
     12 1607  428 (S37W56 L= 20)   484     44        0 (1589/ 397-1627/ 455) -78.7 W53-W57 S39-S34
  

    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