Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1517 1747 (N38W59 L=237)   552     55        1 (1479/1727-1552/1764)   0.0 W54-W65 N36-N38
      2 1571 1485 (N19W45 L=223)   495    165        0 (1510/1437-1644/1529) -48.8 W41-W49 N15-N23
      3  536 1487 (N36E28 L=150)   374     72        0 ( 495/1471- 580/1502)  -8.1 E31-E24 N35-N37
      4  426 1444 (N35E38 L=140)   612     97        0 ( 393/1402- 461/1485)  45.0 E40-E34 N32-N37
      5 1429  986 (S08W24 L=202)   505    183        0 (1377/ 918-1498/1021)  18.4 W19-W28 S11-S05
      6 1007  684 (S20E06 L=172)   646    263        0 ( 915/ 607-1088/ 789)  31.6 E13-W00 S23-S14
      7  234  608 (S13E68 L=110)  2528    185        6 ( 179/ 535- 274/ 678) -58.0 E77-E60 S17-S08
      8  805  621 (S20E20 L=158)   220     54        0 ( 775/ 605- 839/ 634) -33.7 E21-E18 S21-S19
      9 1260  574 (S31W08 L=186)   241     61        0 (1228/ 554-1287/ 602)  26.6 W06-W10 S31-S29
     10  536  462 (S26E45 L=133)   376     86        0 ( 515/ 424- 563/ 507)  71.6 E46-E43 S29-S23
  

    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