Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1029 1422 (N27W12 L=312)   221     57        0 (1002/1400-1059/1444)   0.0 W09-W14 N25-N27
      2 1654 1395 (N07W50 L=350)   364     98        0 (1608/1371-1695/1423) -33.7 W47-W53 N05-N09
      3 1899 1280 (S07W71 L= 11)  1841    145        1 (1870/1221-1919/1363)  18.4 W63-W78 S10-S03
      4  169 1305 (N43E67 L=233)   525     57        0 ( 142/1291- 198/1320)   0.0 E73-E61 N41-N43
      5 1118 1151 (N10W08 L=308)  1331    505        0 ( 918/1025-1305/1239) -46.7 E01-W15 S01-N17
      6 1572 1175 (S02W36 L=335)   627    192        1 (1507/1125-1625/1212) -14.0 W30-W38 S04-N00
      7  594 1000 (N15E25 L=274)   891    194        1 ( 558/ 916- 649/1093)  84.8 E27-E23 N09-N20
      8 1573  748 (S26W26 L=326)  1546    345        0 (1496/ 613-1691/ 876)  25.3 W17-W36 S31-S21
      9  502  801 (N06E36 L=263)   705    119        3 ( 449/ 774- 567/ 838) -29.7 E40-E32 N03-N08
     10 1082  532 (S24E10 L=289)  1627    366        4 ( 966/ 418-1195/ 672)  23.2 E21-W00 S28-S18
     11  538  513 (S11E44 L=255)  1384    305        1 ( 467/ 461- 607/ 600) -74.7 E48-E39 S15-S04
     12 1417  492 (S37W09 L=309)   334     84        0 (1380/ 467-1453/ 526)   9.5 W06-W12 S38-S35
     13  398  464 (S11E59 L=241)   558     64        1 ( 374/ 427- 423/ 494)  78.7 E60-E56 S14-S09
     14  717  445 (S20E36 L=264)   337    114        3 ( 687/ 417- 759/ 490) -68.2 E37-E33 S21-S16
     15  546  355 (S21E54 L=246)  1335    177        3 ( 491/ 298- 615/ 410) -78.7 E58-E48 S26-S16
     16 1277  315 (S46E07 L=293)   359     76        0 (1256/ 281-1296/ 354)  39.8 E09-E03 S48-S42
     17  840  250 (S36E40 L=260)   801     93        1 ( 804/ 215- 880/ 282)  12.5 E44-E35 S38-S33
  

    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