Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  926 1718 (N33E30 L=302)   912    162        1 ( 850/1682- 993/1745)  36.9 E33-E25 N28-N35
      2 1499 1562 (N42W15 L=347)   383     60        0 (1467/1542-1527/1583)  45.0 W12-W17 N39-N44
      3 1307 1097 (N06W14 L=346)   959    319        1 (1226/1015-1441/1187)  61.9 W10-W19 N00-N15
      4  290 1111 (S19E51 L=281)   329     32        0 ( 269/1088- 307/1131) -26.6 E53-E49 S20-S18
      5  592 1097 (S13E27 L=305)   540    181        0 ( 530/1063- 665/1127)  18.4 E30-E21 S16-S11
      6 1730 1009 (N15W44 L= 16)   758    193        1 (1672/ 949-1797/1091) -11.3 W37-W52 N13-N17
      7  252  859 (S36E53 L=279)   504     81        1 ( 222/ 837- 298/ 879)  16.7 E57-E47 S37-S34
      8 1141  806 (S14W13 L=345)   276    117        0 (1092/ 779-1187/ 844)  56.3 W10-W14 S17-S11
      9  616  724 (S34E17 L=315)   526    130        4 ( 561/ 705- 697/ 748)  21.8 E21-E11 S36-S31
     10  411  676 (S44E35 L=297)   723     84        0 ( 366/ 654- 467/ 701)  15.3 E39-E28 S45-S41
     11 1383  451 (S26W43 L= 15)  1974    421        2 (1212/ 403-1580/ 515)  38.7 W33-W54 S34-S15
     12 1150  396 (S37W32 L=  4)   308     59        1 (1120/ 374-1174/ 416)   0.0 W29-W34 S38-S35
  

    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