Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  220 1309 (N35E65 L=151)   501     83        0 ( 180/1287- 251/1322)  -4.1 E73-E59 N33-N35
      2  679 1155 (N20E20 L=196)   621    220        0 ( 591/1126- 772/1188) -18.4 E25-E13 N17-N23
      3 1603  963 (S08W35 L=251)   536    119        1 (1560/ 914-1642/1017) -66.8 W32-W38 S11-S04
      4  804  703 (S09E19 L=197)   475     72        2 ( 776/ 658- 825/ 747)  68.2 E21-E17 S11-S06
      5  318  603 (S10E59 L=156)   842    129        0 ( 288/ 556- 355/ 666) -82.9 E61-E57 S13-S05
      6  366  628 (S08E53 L=162)   502     72        1 ( 337/ 597- 397/ 649)  26.6 E56-E50 S10-S07
      7 1689  572 (S36W43 L=259)  1537    195        2 (1599/ 535-1760/ 621)   0.0 W33-W51 S38-S34
      8  525  510 (S18E44 L=172)   609    189        0 ( 467/ 480- 586/ 552)   0.0 E47-E37 S20-S14
      9 1429  515 (S34W19 L=234)   601    119        1 (1374/ 481-1476/ 551)  14.0 W14-W22 S35-S31
     10 1043  452 (S30E10 L=206)   931    194        0 ( 984/ 382-1117/ 523)  30.3 E15-E03 S34-S25
     11  721  487 (S22E30 L=185)   397     67        1 ( 694/ 466- 760/ 512)  11.3 E32-E27 S23-S20
     12  572  335 (S31E51 L=165)  2037    237        2 ( 538/ 241- 616/ 427)  28.2 E69-E41 S39-S24
  

    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