Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  673 1153 (N20E20 L=196)   573    209        0 ( 592/1124- 765/1183) -14.0 E25-E13 N17-N22
      2 1595  960 (S08W35 L=251)   521    124        1 (1553/ 917-1637/1015) -66.8 W32-W37 S11-S04
      3  797  701 (S09E19 L=197)   442     76        1 ( 770/ 660- 816/ 746)  68.2 E21-E17 S11-S06
      4  334  611 (S09E57 L=159)  1299    185        2 ( 285/ 555- 392/ 662) -76.0 E62-E50 S13-S05
      5 1684  571 (S36W43 L=259)  1522    190        2 (1597/ 534-1753/ 620)   0.0 W33-W50 S38-S33
      6  512  504 (S18E46 L=171)   552    141        0 ( 463/ 480- 564/ 550) -33.7 E48-E40 S19-S14
      7 1419  515 (S33W18 L=234)   681    121        0 (1368/ 481-1469/ 550)  23.2 W13-W22 S35-S31
      8 1035  449 (S30E10 L=206)   991    223        1 ( 961/ 380-1114/ 521)  26.6 E17-E03 S34-S26
      9  716  486 (S22E31 L=185)   462     87        1 ( 682/ 461- 765/ 512)   8.1 E33-E26 S23-S19
     10  567  332 (S31E52 L=164)  2123    231        3 ( 532/ 237- 617/ 418)  27.3 E71-E41 S40-S25
  

    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