Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  655 1700 (N62E08 L=245)   422     84        0 ( 633/1662- 686/1732)  35.5 E10-E02 N58-N63
      2 1271 1637 (N35W38 L=290)  1366    234        2 (1165/1580-1374/1682) -56.3 W33-W42 N28-N40
      3  438 1616 (N62E41 L=212)   465     57        1 ( 411/1597- 461/1640)  23.2 E45-E36 N59-N63
      4  311 1153 (N31E48 L=205)   799    152        1 ( 239/1131- 393/1181) -16.4 E56-E39 N28-N33
      5 1786  892 (S23W47 L=300)  2067    307        3 (1705/ 778-1874/1003)  15.9 W38-W59 S28-S18
      6 1691  721 (S30W36 L=289)   550    105        2 (1644/ 699-1746/ 742) -10.3 W31-W42 S31-S28
      7  887  429 (S25E26 L=227)   700    219        0 ( 824/ 391- 937/ 486) -33.7 E29-E21 S27-S21
      8 1538  420 (S47W20 L=273)   522    106        0 (1510/ 375-1579/ 459)  35.5 W17-W24 S51-S44
      9  556  435 (S17E48 L=205)   417     94        0 ( 505/ 415- 599/ 453) -23.2 E51-E44 S18-S14
     10  519  341 (S23E59 L=194)  3789    217        5 ( 441/ 276- 618/ 414) -45.0 E69-E50 S27-S16
     11  624  382 (S22E46 L=206)   309     70        0 ( 589/ 368- 656/ 396)   0.0 E49-E43 S23-S20
  

    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