Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  973 1881 (N56E40 L= 36)  1596    179        4 ( 877/1851-1052/1910)  33.7 E47-E33 N50-N62
      2  687 1562 (N23E35 L= 40)   413     96        1 ( 643/1539- 736/1583)   0.0 E39-E31 N22-N24
      3  580 1517 (N18E41 L= 35)   424    109        2 ( 554/1475- 610/1568) -35.5 E45-E38 N16-N21
      4 1745 1402 (N38W49 L=125)   455     60        0 (1728/1374-1762/1435) -90.0 W47-W51 N36-N40
      5 1289 1240 (N16W11 L= 87)  1267    404        3 (1138/1102-1396/1320) -19.3 W00-W20 N10-N19
      6 1514 1244 (N22W27 L=102)   363     81        0 (1488/1215-1533/1284)  78.7 W24-W28 N19-N24
      7 1535 1169 (N18W29 L=105)   383    106        2 (1488/1138-1582/1188)  23.2 W25-W32 N15-N19
      8  877  748 (S21E03 L= 73)   292     96        1 ( 835/ 718- 918/ 777)  59.0 E05-E01 S23-S18
      9 1220  553 (S26W24 L=100)  1277    339        2 (1082/ 514-1360/ 644)  47.1 W15-W30 S31-S16
     10 1622  599 (S14W51 L=127)   389     31        1 (1598/ 580-1647/ 620)  -9.5 W48-W54 S14-S12
     11  962  310 (S50W17 L= 93)   307     30        0 ( 942/ 289- 982/ 332) -18.4 W14-W20 S51-S48
     12 1241  241 (S47W47 L=122)  2470    245        4 (1137/ 213-1357/ 296)   9.5 W34-W58 S50-S42
  

    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