Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  993 1649 (N40E14 L=123)   498    133        3 ( 959/1601-1024/1696)  90.0 E17-E12 N35-N44
      2 1403 1568 (N41W20 L=157)   526    158        0 (1352/1524-1453/1602)  37.9 W15-W24 N37-N44
      3  373 1453 (N16E54 L= 84)   536    145        2 ( 323/1417- 418/1503)  54.5 E57-E50 N13-N20
      4 1296 1465 (N32W12 L=149)   221     26        0 (1278/1443-1316/1484)  71.6 W10-W13 N30-N33
      5  938 1259 (N13E09 L=129)   262     66        1 ( 902/1240- 966/1287)  56.3 E10-E07 N11-N14
      6  950 1164 (N07E06 L=131)   688    292        1 ( 845/1053-1054/1232) -26.6 E13-W01 N02-N10
      7 1877 1084 (N16W63 L=201)   607    101        1 (1853/1039-1895/1129) -45.0 W59-W66 N13-N18
      8  618  966 (S09E24 L=113)   403    165        0 ( 544/ 937- 681/1004)  24.0 E29-E20 S11-S06
      9  955  892 (S09E02 L=135)   523     89        4 ( 929/ 829- 986/ 954) -82.9 E04-W00 S12-S04
     10 1557  748 (S09W39 L=177)   370    140        0 (1519/ 687-1586/ 823)  78.7 W37-W40 S13-S03
  

    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