Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  757 1531 (N21E30 L=230)   629    172        0 ( 739/1461- 778/1612) -60.9 E33-E27 N16-N25
      2 1212 1574 (N35E02 L=257)   242     34        1 (1189/1554-1234/1596) -14.0 E04-W00 N34-N36
      3 1064 1482 (N25E09 L=250)   440     86        0 (1028/1445-1102/1522) -23.2 E12-E05 N23-N26
      4  545 1475 (N12E42 L=217)   440     80        1 ( 520/1444- 566/1513) -45.0 E44-E39 N10-N14
      5  625 1378 (N08E33 L=227)   371     63        0 ( 609/1338- 647/1408) -76.0 E34-E31 N06-N10
      6  205 1168 (S12E62 L=197)  1090    114        1 ( 174/1115- 235/1228)  74.1 E65-E59 S15-S08
      7 1310 1187 (N13W13 L=272)   250     79        0 (1267/1172-1351/1202)  31.0 W10-W15 N11-N14
      8 1054 1132 (N04E01 L=259)   383    155        1 ( 990/1102-1101/1163)  51.3 E03-W01 N00-N06
      9  550 1042 (S13E29 L=230)   419    114        2 ( 496/1009- 603/1072)  23.2 E32-E25 S14-S09
     10  172  978 (S24E61 L=198)  1480    235        2 ( 141/ 915- 192/1042) -69.4 E64-E58 S28-S20
     11 1112  950 (S05W07 L=266)   198     44        1 (1086/ 934-1143/ 968)  18.4 W05-W08 S05-S04
     12 1176  761 (S14W15 L=275)   402     73        2 (1140/ 739-1209/ 789) -11.3 W12-W17 S16-S12
     13 1382  719 (S12W29 L=288)   215     70        0 (1356/ 691-1398/ 754)  90.0 W28-W29 S13-S09
     14  536  702 (S34E25 L=234)   647    132        1 ( 508/ 643- 574/ 745)  69.4 E27-E22 S38-S30
     15  324  701 (S39E44 L=215)   705    148        0 ( 293/ 658- 363/ 741)  45.0 E47-E40 S42-S35
     16  335  597 (S47E45 L=215)   706     84        0 ( 307/ 556- 368/ 643) -56.3 E47-E41 S49-S43
     17 1089  551 (S29W17 L=276)   401     75        1 (1041/ 533-1131/ 574)  18.4 W13-W19 S31-S27
     18 1619  480 (S19W57 L=316)   413     43        0 (1590/ 465-1645/ 498)  31.0 W54-W59 S20-S17
     19  855  458 (S42W05 L=264)   288     44        1 ( 832/ 440- 882/ 484)  33.7 W03-W06 S43-S39
  

    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