Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  915 1762 (N39E34 L=297)   858    138        0 ( 888/1700- 945/1827) -34.7 E41-E27 N34-N43
      2 1465 1715 (N56W13 L=344)  1200    110        2 (1412/1688-1524/1740)  34.7 W06-W20 N50-N60
      3  637 1676 (N25E47 L=284)   667    129        2 ( 605/1623- 671/1739) -24.8 E54-E41 N21-N28
      4 1266 1273 (N16W07 L=338)   275     83        0 (1227/1249-1307/1290)   9.5 W04-W10 N15-N17
      5  167 1009 (S25E66 L=265)  1522    216        3 ( 117/ 918- 241/1070)  42.5 E77-E55 S31-S20
      6  493  764 (S33E28 L=303)   690    190        2 ( 408/ 736- 595/ 786)  23.6 E36-E20 S35-S28
      7 1007  620 (S27W10 L=341)   590    132        0 ( 978/ 557-1043/ 678) -83.7 W09-W11 S32-S23
      8 1360  482 (S25W39 L= 10)   383    118        0 (1293/ 456-1403/ 509)  45.0 W35-W41 S28-S22
      9  402  436 (S61E45 L=286)  1049    108        2 ( 363/ 401- 440/ 483) -90.0 E54-E37 S63-S56
     10  502  383 (S62E29 L=302)   754     59        1 ( 476/ 353- 525/ 417) -80.5 E33-E24 S65-S58
     11  664  338 (S60E05 L=326)   376     29        0 ( 637/ 320- 685/ 355)   0.0 E08-E02 S61-S57
     12  868  234 (S60W24 L=355)   447     52        0 ( 849/ 209- 893/ 267) -21.8 W18-W28 S62-S57
  

    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