Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1759 1487 (N31W66 L=337)   971     90        2 (1710/1463-1802/1513) -17.1 W60-W74 N29-N33
      2  554 1462 (N31E36 L=235)   332    102        1 ( 513/1439- 607/1489) -12.5 E39-E30 N29-N32
      3  966 1450 (N30E04 L=267)   303    143        0 ( 907/1419-1031/1488) -26.6 E08-W00 N28-N33
      4  143 1206 (N12E75 L=196)  1661    124        0 ( 120/1165- 166/1264)  63.4 E79-E69 N09-N15
      5 1836  964 (S02W60 L=331)  4041    410       13 (1758/ 835-1916/1122) -64.8 W52-W72 S09-N08
      6  190  938 (S04E63 L=207)  2539    260        8 ( 136/ 875- 269/1033) -19.7 E72-E54 S08-N01
      7  611  987 (N01E26 L=245)   304    133        0 ( 575/ 932- 662/1031) -39.8 E28-E22 S02-N03
      8  256  822 (S11E57 L=214)   934    142        1 ( 231/ 750- 275/ 883)  66.0 E60-E53 S15-S06
      9 1743  666 (S20W55 L=326)   468     84        3 (1718/ 639-1760/ 698)  76.0 W52-W57 S22-S18
     10  331  634 (S23E53 L=217)  1045    127        2 ( 296/ 577- 382/ 673)  26.6 E59-E46 S26-S19
     11  402  586 (S26E48 L=223)   550    133        0 ( 357/ 552- 443/ 634)  22.6 E53-E41 S28-S22
     12 1632  561 (S27W47 L=318)   664    148        0 (1607/ 491-1656/ 627) -68.2 W43-W49 S32-S22
     13 1366  521 (S29W25 L=296)  1751    319        6 (1232/ 459-1500/ 591)  20.2 W15-W34 S33-S24
  

    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