Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1253 1852 (N65W43 L=349)  1240    121        1 (1188/1833-1308/1884)   5.2 W31-W65 N62-N67
      2 1543 1574 (N37W47 L=353)   402     63        0 (1507/1559-1576/1591)   7.1 W43-W51 N36-N38
      3  942 1380 (N27E04 L=301)  1056    384        2 ( 882/1290-1048/1471) -57.5 E08-W02 N20-N33
      4 1544 1384 (N24W39 L=345)   275     65        0 (1520/1354-1566/1414) -63.4 W37-W40 N22-N26
      5 1461 1225 (N15W30 L=336)   262     57        0 (1439/1191-1485/1250)  45.0 W27-W31 N12-N16
      6 1175 1139 (N11W10 L=316)   191     33        1 (1151/1121-1200/1156) -18.4 W08-W11 N09-N11
      7 1746  857 (S10W51 L=357)   595     33        0 (1719/ 831-1772/ 880)  36.9 W48-W53 S11-S08
      8  782  813 (S08E16 L=290)   754    194        0 ( 708/ 772- 844/ 849)  -7.1 E20-E12 S10-S05
      9  328  650 (S19E54 L=251)   526     88        0 ( 301/ 617- 366/ 674)  20.6 E57-E49 S21-S17
     10 1648  609 (S26W46 L=352)   409     66        1 (1623/ 581-1677/ 644) -29.7 W42-W49 S27-S23
     11 1539  592 (S26W36 L=342)   354     92        0 (1512/ 572-1577/ 629) -33.7 W33-W39 S27-S23
     12 1171  149 (S68W17 L=323)   664     47        0 (1159/ 118-1181/ 179) -47.3 W13-W27 S75-S62
  

    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