Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  920 1201 (N18E01 L=325)   362    120        0 ( 862/1182- 987/1229) -23.2 E04-W03 N16-N20
      2 1793 1030 (S17W49 L= 15)   405     90        2 (1772/ 991-1819/1066)  38.7 W46-W52 S18-S14
      3  133 1020 (N26E64 L=262)   701    121        1 ( 110/ 978- 156/1056) -68.2 E66-E61 N23-N28
      4  675  886 (N07E23 L=303)   849    230        1 ( 565/ 848- 776/ 930) -30.3 E29-E17 N03-N11
      5  317  854 (N13E50 L=277)  1696    546        5 ( 174/ 779- 456/ 914)  -4.1 E65-E37 N09-N15
      6 1350  726 (S20W11 L=337)   272     92        0 (1332/ 678-1365/ 772)  51.3 W08-W12 S22-S17
      7  858  427 (S24E28 L=298)  1591    413        3 ( 819/ 288- 894/ 578)  45.0 E37-E20 S33-S15
      8  315  450 (S11E71 L=256)  4570    298        4 ( 255/ 363- 396/ 526) -71.6 E79-E66 S18-S05
      9  422  394 (S16E62 L=265)  1393    241        6 ( 373/ 348- 477/ 453) -63.4 E65-E57 S20-S12
     10 1189  333 (S41E12 L=314)   269     79        1 (1167/ 305-1219/ 359) -90.0 E13-E11 S44-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