Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  732 1662 (N41E31 L=289)  4083    542        7 ( 489/1539- 937/1772) -10.8 E74-E11 N33-N47
      2 1181 1547 (N37W07 L=327)  1202    238        2 (1060/1517-1282/1580)  10.6 E01-W15 N33-N39
      3 1371 1466 (N32W22 L=341)   414    135        0 (1325/1411-1421/1514) -45.0 W18-W25 N28-N35
      4  250 1425 (N20E67 L=252)  1823    185        7 ( 205/1362- 295/1504) -78.7 E72-E63 N15-N25
      5  439 1445 (N23E47 L=273)   312     47        0 ( 420/1427- 463/1468) -26.6 E48-E44 N21-N24
      6 1235 1346 (N23W11 L=331)   477    166        1 (1208/1275-1250/1427)  78.7 W09-W12 N18-N28
      7 1737 1365 (N28W53 L= 13)   337     56        0 (1718/1334-1750/1389) -90.0 W51-W55 N25-N29
      8  358 1146 (N04E46 L=274)  2081    316        3 ( 329/ 997- 392/1294) -83.7 E48-E43 S05-N13
      9 1848 1241 (N20W64 L= 23)  1145     86        2 (1816/1202-1885/1279) -39.8 W59-W69 N17-N22
     10 1897 1153 (N15W70 L= 29)   663     72        1 (1866/1129-1930/1173)  -4.8 W64-W76 N13-N15
     11  188  971 (S08E62 L=257)   531    123        0 ( 151/ 920- 221/1019)  35.0 E68-E58 S11-S04
     12 1613  916 (S01W39 L=359)   251     56        0 (1592/ 900-1643/ 938)   0.0 W37-W41 S01-N00
     13  573  769 (S17E28 L=292)   274    122        0 ( 521/ 748- 620/ 805) -23.2 E31-E24 S18-S15
     14  377  676 (S25E45 L=275)   367     78        0 ( 349/ 651- 411/ 701)  23.2 E48-E41 S27-S23
     15 1336  650 (S19W23 L=343)  1016    205        1 (1239/ 617-1433/ 691)  22.6 W17-W29 S21-S16
     16 1039  517 (S30W05 L=325)   691    148        1 ( 975/ 498-1112/ 538)  12.5 W00-W09 S31-S28
     17  675  373 (S45E24 L=295)   294     26        1 ( 656/ 356- 694/ 394)  56.3 E26-E22 S46-S43
  

    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