Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  739 1764 (N56E11 L=274)   421     74        0 ( 706/1734- 770/1784)  15.9 E14-E07 N53-N57
      2 1865  962 (S16W61 L=346)   853     60        2 (1835/ 932-1894/ 990) -26.6 W57-W65 S17-S13
      3 1750  923 (S16W48 L=333)   289     41        2 (1729/ 900-1773/ 949) -71.6 W45-W50 S17-S14
      4 1825  871 (S21W56 L=341)   382     76        0 (1806/ 844-1854/ 904) -33.7 W53-W60 S22-S18
      5  668  538 (S23E32 L=253)   242     36        1 ( 644/ 514- 691/ 563) -90.0 E33-E30 S25-S21
      6 1656  502 (S43W45 L=329)   631     87        1 (1629/ 473-1692/ 543) -78.7 W41-W49 S45-S40
      7  665  469 (S28E35 L=250)   402     96        1 ( 633/ 417- 690/ 506)  80.5 E36-E32 S31-S25
      8  964  463 (S33E14 L=271)   494     99        1 ( 931/ 428- 989/ 505)  45.0 E16-E11 S35-S30
      9 1272  437 (S41W08 L=293)   240     27        0 (1255/ 417-1285/ 465)  63.4 W06-W09 S42-S38
     10 1108  316 (S47E08 L=277)   660     95        0 (1053/ 295-1158/ 338) -32.0 E11-E03 S49-S44
  

    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