Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  596 1702 (N61E14 L=123)   655    158        0 ( 527/1670- 657/1747)   9.5 E24-E05 N59-N63
      2  534 1500 (N46E22 L=115)   344     53        1 ( 511/1469- 559/1523)  56.3 E24-E19 N44-N48
      3  960 1482 (N33W09 L=147)   281     86        0 ( 918/1456- 993/1506) -36.9 W06-W11 N30-N34
      4 1182 1410 (N22W20 L=158)   277    121        1 (1136/1377-1235/1434)   0.0 W16-W24 N21-N23
      5  435 1122 (N24E35 L=103)   239     89        0 ( 425/1081- 452/1167)  80.5 E35-E34 N21-N27
      6 1806  966 (S20W50 L=188)   632    154        0 (1766/ 901-1862/1013)  14.0 W45-W57 S22-S18
      7 1759  772 (S31W44 L=181)  1131    312        1 (1652/ 680-1846/ 886)  20.2 W32-W54 S33-S25
      8  914  669 (S12E15 L=122)   214     56        1 ( 894/ 639- 928/ 703)  63.4 E16-E13 S14-S10
      9  566  589 (S09E39 L= 98)   641    106        3 ( 536/ 532- 602/ 644)  35.5 E43-E35 S11-S05
     10 1376  468 (S38W07 L=145)   305     60        2 (1348/ 452-1408/ 490)   0.0 W04-W09 S39-S37
     11  631  341 (S26E48 L= 89)  1750    195        3 ( 582/ 287- 722/ 418) -66.0 E54-E42 S30-S20
     12 1021  360 (S35E21 L=117)   400    104        0 ( 982/ 330-1058/ 403)  14.0 E24-E16 S36-S32
  

    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