Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  309 1536 (N35E68 L=193)   916     90        2 ( 295/1494- 326/1593) -26.6 E75-E62 N32-N38
      2 1103 1447 (N30W06 L=267)   459    175        0 (1022/1416-1180/1484) -20.0 W00-W11 N28-N32
      3  198 1203 (N13E65 L=196)  1621    175        4 ( 156/1140- 228/1289)  63.4 E69-E61 N08-N18
      4 1889 1027 (N01W68 L=329)  2974    218       11 (1819/ 920-1943/1123) -56.3 W58-W79 S05-N07
      5  746  997 (N02E17 L=244)   246    108        1 ( 708/ 957- 789/1025) -31.0 E19-E14 N00-N03
      6  234  928 (S04E58 L=203)   684    127        0 ( 211/ 871- 266/ 978)  49.4 E61-E54 S07-N00
      7  360  803 (S11E46 L=215)   701    176        1 ( 309/ 737- 390/ 872)  66.0 E52-E43 S15-S06
      8  179  736 (S17E71 L=190)   853    101        1 ( 157/ 689- 206/ 776)  22.4 E81-E64 S20-S13
      9  444  617 (S23E43 L=218)  1453    228        4 ( 390/ 529- 507/ 695)  45.0 E50-E36 S29-S17
     10 1721  559 (S28W57 L=318)   681    123        0 (1693/ 500-1741/ 616) -69.4 W54-W59 S32-S24
     11 1466  517 (S30W33 L=294)  2221    449        5 (1232/ 449-1615/ 604)  13.6 W15-W44 S34-S24
     12  439  402 (S40E55 L=206)  1186     87        2 ( 405/ 374- 481/ 432)  23.2 E62-E48 S42-S36
  

    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