Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1165 1761 (N44E18 L=197)   634    192        0 (1103/1682-1216/1825) -23.6 E26-E10 N39-N47
      2  320 1447 (N00E60 L=155)  1065     86        2 ( 296/1400- 350/1485) -76.0 E64-E56 S02-N02
      3 1170 1380 (N18E02 L=213)   779    192        0 (1099/1304-1237/1434) -16.7 E07-W03 N14-N19
      4  357 1347 (S05E52 L=163)   514     83        0 ( 333/1307- 381/1381)  78.7 E54-E49 S07-S02
      5  207 1071 (S24E60 L=155)   475     46        1 ( 188/1042- 224/1104) -76.0 E62-E57 S25-S21
      6  763  916 (S19E12 L=203)  1086    430        0 ( 567/ 820- 930/ 973)  41.8 E22-E03 S29-S12
      7  391  813 (S35E37 L=177)   714    176        1 ( 359/ 738- 451/ 871)  59.0 E40-E32 S40-S30
      8  887  758 (S25E00 L=215)   566    177        0 ( 822/ 719- 972/ 786)  53.1 E03-W04 S28-S20
      9  961  632 (S30W09 L=223)   901    294        1 ( 820/ 613-1085/ 660)  31.0 W00-W15 S34-S25
     10 1264  503 (S28W32 L=247)   491    129        1 (1205/ 484-1322/ 533)   6.3 W27-W36 S28-S25
     11 1555  463 (S20W55 L=270)  1014    144        2 (1485/ 423-1623/ 508)  56.3 W51-W59 S24-S15
     12  377  460 (S62E47 L=168)  1064    142        2 ( 316/ 431- 429/ 494)  -2.7 E60-E36 S64-S59
  

    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