Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  389 1509 (N34E54 L=195)   412     55        2 ( 365/1479- 404/1541) -35.5 E57-E50 N31-N36
      2  812 1455 (N31E15 L=234)   481    136        0 ( 755/1430- 895/1485) -15.3 E19-E08 N29-N33
      3  299 1230 (N15E53 L=195)  1513    311        1 ( 265/1122- 321/1377) -79.4 E58-E50 N08-N24
      4 1572 1129 (N09W36 L=285)  1401    516        7 (1456/ 948-1646/1370) -68.2 W30-W41 S02-N24
      5  944  996 (N02E05 L=244)   314    120        0 ( 886/ 962- 985/1024) -26.6 E08-E02 N00-N03
      6  156  895 (S06E69 L=180)  1337    143        2 ( 129/ 830- 187/ 959) -49.4 E73-E65 S09-S02
      7  510  764 (S13E35 L=214)   552    148        1 ( 459/ 722- 549/ 821)  45.0 E38-E31 S15-S09
      8  246  651 (S21E64 L=185)  1990    183        2 ( 209/ 570- 291/ 734)  40.6 E72-E54 S27-S15
      9  811  690 (S17E14 L=235)   456    158        1 ( 749/ 662- 872/ 706)   7.1 E18-E10 S18-S16
     10  606  595 (S24E30 L=219)  1201    293        2 ( 560/ 500- 665/ 678)  47.5 E35-E24 S30-S18
     11 1677  558 (S28W52 L=300)  1735    203        2 (1574/ 486-1742/ 624)  36.0 W43-W58 S32-S23
     12  536  369 (S42E45 L=204)  1896    132        3 ( 487/ 317- 602/ 420)  27.8 E57-E37 S46-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