Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1117 1789 (N47W38 L=209)   309     49        0 (1100/1770-1134/1813)  21.8 W35-W40 N45-N48
      2 1607 1624 (N20W59 L=231)  2068    301        2 (1499/1581-1722/1656) -34.5 W51-W68 N14-N25
      3  397 1146 (N28E37 L=134)   961    158        1 ( 352/1102- 459/1198) -58.0 E40-E32 N23-N31
      4  165  905 (N18E64 L=107)   607    118        2 ( 143/ 879- 192/ 943)  63.4 E66-E60 N16-N20
      5 1143  564 (S24E06 L=165)   667    129        1 (1098/ 511-1190/ 616)  26.6 E10-E02 S26-S21
      6  429  544 (S08E52 L=119)   352     53        0 ( 407/ 523- 453/ 566)  11.3 E54-E49 S08-S06
      7  621  506 (S14E39 L=132)   255     61        0 ( 585/ 494- 660/ 517) -26.6 E41-E37 S15-S12
      8 1504  444 (S45W16 L=187)   421     45        0 (1474/ 428-1540/ 460) -15.9 W12-W19 S47-S43
  

    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