Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1306 1784 (N51E08 L=238)   665     77        0 (1263/1758-1348/1807)   6.3 E12-E03 N49-N53
      2  405 1617 (N13E62 L=185)   608     77        0 ( 372/1592- 440/1646)  59.0 E64-E59 N10-N15
      3 1677 1434 (N38W33 L=279)   425     46        1 (1656/1412-1702/1458) -26.6 W30-W35 N36-N40
      4  281 1428 (S01E62 L=184)   565     77        3 ( 253/1386- 301/1456) -23.2 E66-E58 S03-N00
      5  232 1329 (S08E63 L=183)   555     77        0 ( 216/1288- 248/1370) -68.2 E65-E61 S10-S05
      6  653 1164 (S08E25 L=222)   562    200        0 ( 621/1082- 721/1227)  84.3 E26-E22 S13-S03
      7 1339 1183 (N11W13 L=259)   448    106        1 (1289/1147-1380/1226) -15.9 W09-W16 N09-N12
      8  142 1168 (S18E70 L=176)  1711    123        2 (  92/1130- 209/1195)  11.3 E81-E61 S20-S15
      9 1001 1087 (S03E03 L=243)   202     80        0 ( 977/1065-1017/1132) -45.0 E04-E01 S04-S01
     10  464  957 (S24E33 L=213)   732    197        1 ( 414/ 874- 541/1021)  59.0 E36-E28 S29-S18
     11 1281  779 (S12W21 L=268)   343    132        1 (1220/ 755-1351/ 805)  39.8 W18-W24 S14-S09
     12 1398  659 (S15W33 L=279)   285     73        0 (1367/ 632-1416/ 693)  90.0 W31-W34 S16-S12
     13  540  538 (S49E20 L=226)  1225    272        2 ( 418/ 508- 662/ 598)   6.8 E33-E08 S51-S44
     14 1663  562 (S11W56 L=302)   330     79        1 (1649/ 519-1688/ 591) -63.4 W53-W58 S13-S09
     15 1218  518 (S28W28 L=274)   406    111        0 (1169/ 500-1271/ 546)  15.9 W24-W31 S29-S26
     16 1582  406 (S22W61 L=307)   965    154        2 (1525/ 356-1640/ 445)  66.8 W57-W64 S25-S18
     17 1471  407 (S26W52 L=298)   528     60        0 (1439/ 387-1511/ 429)   7.1 W48-W56 S27-S24
     18 1076  274 (S48W37 L=283)  1280    267        0 ( 955/ 216-1182/ 335)   2.1 W23-W50 S49-S46
  

    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