Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1460 1628 (N29W45 L=358)  2557    420        3 (1259/1558-1614/1688) -32.5 W32-W57 N21-N36
      2 1724 1505 (N16W62 L= 14)   538     33        0 (1697/1489-1753/1522)   0.0 W58-W65 N14-N16
      3 1406 1457 (N19W33 L=345)   377     86        1 (1357/1437-1449/1476) -38.7 W30-W35 N17-N21
      4  768 1318 (N23E10 L=302)   211     36        0 ( 751/1292- 789/1343)  56.3 E11-E08 N21-N24
      5  486 1263 (N25E31 L=282)  2097    419        2 ( 305/1210- 669/1328)  -5.7 E46-E16 N22-N29
      6 1838 1270 (S00W62 L= 15)  1045    129        0 (1788/1227-1868/1328)  12.1 W55-W69 S01-N02
      7 1626 1098 (S06W39 L=351)   280     36        1 (1597/1080-1652/1114) -26.6 W36-W40 S06-S04
      8 1861  869 (S24W59 L= 11)  1383    161        1 (1831/ 787-1890/ 946)  84.3 W55-W63 S29-S19
      9 1085  768 (S14E01 L=311)   292     82        0 (1068/ 722-1108/ 810)  63.4 E02-W00 S16-S12
     10 1259  745 (S19W09 L=321)   370    108        0 (1225/ 704-1280/ 799)  45.0 W06-W10 S21-S16
     11 1382  698 (S24W16 L=328)   269     91        0 (1362/ 653-1404/ 745)  59.0 W14-W18 S26-S21
     12  440  435 (S23E56 L=256)   952    179        0 ( 400/ 366- 463/ 521)  50.7 E64-E52 S27-S16
  

    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