Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  366 1593 (N10E64 L= 16)  2311    357        5 ( 300/1477- 421/1678) -77.9 E67-E58 N02-N16
      2  444 1564 (N10E54 L= 26)   363     90        3 ( 419/1519- 470/1597)  90.0 E55-E52 N06-N12
      3 1222 1368 (N18W02 L= 82)   683    267        0 (1203/1252-1255/1492) -61.7 E02-W06 N12-N26
      4 1565 1380 (N30W23 L=103)   278     41        1 (1543/1359-1589/1403) -33.7 W21-W25 N29-N31
      5  217 1311 (S09E65 L= 16)  1758    338        7 ( 164/1251- 273/1364) -54.5 E71-E57 S12-S05
      6  735 1284 (N00E23 L= 57)   358    121        1 ( 678/1247- 801/1306)   0.0 E26-E18 S01-N01
      7 1413 1236 (N16W16 L= 96)   386    137        0 (1393/1165-1430/1302) -54.5 W13-W18 N12-N19
      8  857 1068 (S09E10 L= 70)   201     63        1 ( 831/1043- 894/1084)  33.7 E11-E08 S10-S07
      9 1228  932 (S06W14 L= 94)   456    108        1 (1208/ 880-1247/ 991) -74.1 W12-W14 S09-S02
     10  165  963 (S31E60 L= 21)  1577    427        3 (  88/ 935- 241/ 988)  10.6 E68-E52 S32-S27
     11  653  890 (S24E18 L= 62)   705    179        3 ( 591/ 823- 703/ 950)  68.2 E20-E15 S29-S19
     12 1821  853 (N09W57 L=138)   635    108        3 (1801/ 810-1842/ 899) -51.3 W53-W61 N06-N11
     13  410  763 (S38E34 L= 46)  1350    335        4 ( 308/ 728- 530/ 813)  21.8 E44-E24 S42-S33
     14 1824  692 (N01W66 L=146)  1490    354        4 (1764/ 666-1897/ 748)  22.4 W59-W77 S02-N05
     15  905  621 (S32W05 L= 86)   840    252        0 ( 825/ 542-1017/ 732)  76.8 W01-W08 S39-S22
     16 1253  624 (S22W26 L=106)   540    179        1 (1164/ 572-1327/ 694)  65.6 W21-W27 S26-S15
     17 1527  602 (S14W44 L=124)   305    129        1 (1481/ 576-1564/ 635)   0.0 W39-W47 S15-S12
     18 1269  473 (S29W34 L=114)   272     59        0 (1237/ 460-1311/ 487)  31.0 W31-W36 S30-S27
     19 1424  323 (S32W53 L=134)   426     44        1 (1393/ 308-1453/ 343)  36.9 W51-W55 S33-S30
  

    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