Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1360 1583 (N32W38 L= 30)   501    100        0 (1304/1569-1406/1599) -14.0 W33-W41 N30-N33
      2  719 1466 (N40E11 L=340)   596    118        1 ( 662/1420- 775/1508)  16.7 E16-E06 N37-N42
      3  681 1292 (N30E17 L=335)   219     78        1 ( 650/1279- 721/1307) -26.6 E18-E13 N28-N31
      4 1058 1259 (N20W07 L=359)   749    269        3 ( 914/1229-1183/1282) -25.0 E01-W14 N16-N23
      5 1358 1211 (N11W24 L= 16)   430    132        1 (1285/1187-1412/1225) -20.6 W19-W27 N08-N12
      6 1051  591 (S19E08 L=344)   379    154        0 (1015/ 519-1084/ 644)  45.0 E11-E04 S23-S16
      7  533  536 (S13E43 L=309)   399     86        0 ( 514/ 503- 564/ 578) -78.7 E44-E41 S15-S10
      8  892  323 (S35E27 L=324)   478     74        1 ( 867/ 290- 924/ 343)  31.0 E29-E24 S37-S32
      9 1275  255 (S52E02 L=350)   781    139        2 (1214/ 235-1345/ 273) -20.0 E07-W05 S54-S49
  

    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