Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1410 1748 (N38W52 L= 68)   548     75        1 (1383/1713-1447/1779)  68.2 W48-W55 N34-N39
      2 1571 1719 (N32W65 L= 81)  1503    130        2 (1511/1691-1625/1752) -24.4 W57-W71 N28-N34
      3  612 1700 (N55E20 L=357)  1530    178        3 ( 566/1644- 666/1745) -33.7 E25-E13 N50-N59
      4  916 1600 (N39W06 L= 22)  1412    290        2 ( 801/1535-1054/1647) -45.0 E01-W13 N31-N45
      5 1045 1394 (N23W09 L= 25)   427     64        1 (1009/1369-1081/1417)  11.3 W06-W11 N21-N24
      6  368 1273 (N29E41 L=336)  1521    351        2 ( 267/1181- 461/1385) -38.7 E53-E32 N21-N38
      7 1190 1289 (N14W15 L= 31)   510    104        0 (1126/1260-1238/1319) -39.8 W11-W17 N11-N16
      8  203 1083 (N20E58 L=318)  2582    354        2 (  94/ 981- 324/1210) -22.1 E83-E46 N13-N29
      9 1903 1128 (S10W67 L= 84)  1270    175        0 (1884/1037-1931/1204)  45.0 W62-W73 S14-S05
     10  432 1096 (N17E36 L=340)   853    146        2 ( 382/1058- 482/1133)   7.1 E40-E32 N14-N19
     11  462  807 (S00E39 L=337)   535     93        1 ( 418/ 783- 510/ 835) -18.4 E41-E35 S01-N01
     12 1064  670 (S19E05 L= 12)   200     69        0 (1045/ 640-1078/ 707)  45.0 E06-E03 S21-S17
     13 1309  630 (S27W10 L= 26)   364     73        2 (1288/ 586-1331/ 667)  45.0 W07-W11 S29-S24
     14  520  526 (S18E44 L=332)   317     40        1 ( 498/ 500- 544/ 544) -71.6 E46-E42 S19-S16
     15  325  455 (S19E69 L=307)  1634    108        5 ( 286/ 412- 363/ 496) -33.7 E76-E63 S22-S16
     16 1061  401 (S37E12 L=  4)   289     31        1 (1041/ 377-1084/ 422)  33.7 E13-E09 S38-S35
  

    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