Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1603 1678 (N28W70 L= 21)  1644    152        2 (1542/1620-1658/1722) -90.0 W63-W74 N23-N31
      2  593 1546 (N50E20 L=291)   317     43        0 ( 566/1527- 616/1567)  18.4 E23-E17 N48-N51
      3 1370 1523 (N27W37 L=348)   331    110        0 (1326/1475-1408/1556) -71.6 W35-W38 N23-N29
      4 1506 1523 (N23W47 L=358)   327     67        2 (1483/1486-1532/1549)  26.6 W44-W50 N20-N23
      5 1630 1508 (N18W56 L=  8)  1017    176        1 (1552/1460-1702/1544) -48.4 W51-W60 N13-N22
      6  840 1453 (N37E01 L=310)   280     42        1 ( 819/1422- 860/1476)  76.0 E02-W00 N34-N38
      7 1020 1254 (N20W05 L=317)   218     59        0 (1006/1215-1033/1284)  53.1 W03-W06 N18-N22
      8 1755 1159 (S04W51 L=  3)   387     81        1 (1734/1118-1784/1189)  26.6 W48-W54 S06-S02
      9  193 1078 (N26E63 L=248)   661     64        3 ( 162/1045- 219/1112) -35.5 E68-E59 N23-N28
     10  276 1042 (N23E53 L=259)  1699    370        2 ( 173/ 935- 412/1106)   8.4 E66-E39 N17-N26
     11 1562  990 (S08W32 L=343)   276     42        0 (1544/ 967-1581/1015)  18.4 W29-W33 S09-S06
     12 1845  854 (S26W56 L=  7)  1843    185        2 (1781/ 792-1918/ 917)   6.0 W47-W68 S28-S22
     13 1698  865 (S20W40 L=351)   805    231        1 (1580/ 832-1766/ 907) -14.0 W30-W46 S22-S15
     14 1067  807 (S06E03 L=309)   419    144        2 (1024/ 736-1102/ 860)  45.0 E06-W00 S09-S03
     15 1511  744 (S21W23 L=335)   290    117        0 (1469/ 712-1565/ 771)   8.1 W20-W27 S22-S20
     16  501  707 (S01E40 L=271)   406    115        1 ( 487/ 650- 515/ 766)  56.3 E41-E37 S04-N02
     17 1041  339 (S36E18 L=293)  3349    648        8 ( 891/ 135-1205/ 639)  32.6 E50-W01 S52-S19
     18 1396  458 (S38W11 L=322)   356     39        0 (1376/ 428-1410/ 494)  78.7 W08-W12 S40-S35
     19  572  363 (S24E51 L=260)   383     95        0 ( 524/ 348- 617/ 379)   0.0 E56-E46 S24-S23
  

    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