Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  921 1830 (N64E14 L= 10)   935    136        3 ( 862/1806- 992/1849)  -3.2 E21-E03 N61-N66
      2  679 1482 (N33E26 L=358)   273     37        1 ( 658/1459- 702/1500) -21.8 E28-E23 N31-N34
      3 1372 1398 (N27W25 L= 49)   392    146        1 (1329/1347-1419/1440) -50.2 W22-W27 N23-N29
      4 1019 1327 (N23E00 L= 24)   297    117        0 ( 977/1283-1063/1366) -50.2 E03-W02 N19-N25
      5  321 1289 (N19E53 L=331)   334     61        0 ( 296/1268- 352/1312)  31.0 E55-E50 N17-N20
      6 1828  852 (S09W61 L= 85)   573    120        0 (1807/ 800-1839/ 910)  81.9 W59-W62 S12-S05
      7 1750  855 (S08W52 L= 76)   493     15        1 (1728/ 831-1773/ 877)  45.0 W49-W54 S09-S06
      8 1671  838 (S09W45 L= 69)   382     58        2 (1647/ 802-1696/ 873) -59.0 W42-W46 S11-S06
      9 1610  756 (S14W40 L= 64)   792    217        2 (1559/ 715-1670/ 810)  18.4 W36-W45 S16-S10
     10 1059  696 (S17W02 L= 26)   748    209        3 ( 975/ 656-1153/ 731)  20.0 E03-W08 S19-S14
     11 1541  581 (S26W38 L= 62)  1005    178        0 (1468/ 546-1595/ 647)  40.6 W31-W41 S28-S21
     12  522  488 (S32E40 L=344)  2162    348        4 ( 442/ 415- 621/ 573) -39.8 E45-E32 S37-S25
  

    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