Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  800 1537 (N22E22 L=346)   219     54        0 ( 782/1507- 815/1565) -53.1 E23-E20 N19-N23
      2  485 1518 (N17E44 L=324)  1047    186        5 ( 411/1488- 580/1544)   8.7 E49-E36 N15-N19
      3  867 1470 (N18E16 L=352)   532    102        1 ( 808/1438- 912/1496)   0.0 E19-E12 N16-N19
      4  979 1441 (N18E09 L=359)   405    176        2 ( 944/1407-1022/1496)  63.4 E10-E06 N15-N21
      5  845 1367 (N12E15 L=353)   737    208        3 ( 784/1302- 899/1445)  77.5 E18-E12 N07-N16
      6  675 1361 (N09E26 L=342)   231    112        0 ( 652/1304- 693/1403) -51.3 E27-E23 N06-N11
      7  951 1369 (N13E09 L=359)   276     83        2 ( 926/1335- 992/1390)  26.6 E10-E06 N11-N14
      8  653 1015 (S11E22 L=346)   425    155        1 ( 598/ 968- 692/1074)  66.8 E25-E19 S14-S07
      9 1154 1063 (S02W07 L= 15)   254    133        0 (1085/1041-1214/1073)  23.2 W03-W10 S04-N00
     10 1533  956 (S02W32 L= 40)   384    164        2 (1516/ 885-1550/1038) -82.9 W30-W34 S06-N02
     11  166  893 (S22E64 L=304)  1705    236        5 (  92/ 819- 226/ 963)  20.6 E80-E55 S26-S17
     12  606  782 (S26E24 L=344)  1284    326        2 ( 457/ 739- 757/ 836)  10.3 E35-E13 S29-S21
     13 1339  682 (S21W25 L= 33)  1152    416        2 (1220/ 566-1506/ 807)  57.7 W20-W33 S30-S11
     14 1073  763 (S21W07 L= 15)   435    192        0 ( 984/ 736-1162/ 783)  16.7 W01-W12 S21-S18
     15 1562  633 (S20W42 L= 50)   843    163        4 (1525/ 563-1589/ 718) -66.0 W38-W44 S24-S15
     16 1001  552 (S35W06 L= 14)   915    250        1 ( 880/ 532-1115/ 587)   7.1 E02-W14 S36-S32
     17  396  556 (S43E45 L=323)   768    229        2 ( 298/ 531- 522/ 582)   0.0 E57-E31 S44-S41
     18 1408  285 (S45W52 L= 60)  2108    230        5 (1316/ 237-1495/ 339)  -2.0 W39-W67 S47-S41
  

    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