Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  568 1542 (N12E42 L=189)  1088    299        1 ( 523/1405- 642/1661) -80.5 E44-E39 N03-N21
      2 1325 1236 (N12W12 L=243)  1947    438        1 (1231/1040-1452/1456) -39.3 W01-W23 N04-N23
      3 1488 1418 (N28W18 L=249)   255     49        0 (1464/1391-1505/1442) -45.0 W16-W19 N26-N29
      4  687 1324 (N02E27 L=204)  1473    433        4 ( 486/1231- 811/1368)   2.5 E40-E17 S01-N05
      5  237 1277 (S09E60 L=171)   387    102        0 ( 206/1247- 265/1317)  38.7 E63-E58 S10-S06
      6  788 1132 (S07E16 L=215)   664    183        0 ( 746/1073- 813/1219) -77.5 E18-E13 S10-S01
      7 1793 1004 (N13W49 L=281)  2754    404        4 (1758/ 812-1819/1214) -59.4 W43-W57 N03-N25
      8  407 1139 (S14E41 L=190)   253     69        0 ( 377/1111- 430/1155) -11.3 E43-E38 S15-S13
      9  623 1082 (S13E25 L=206)   866    210        1 ( 532/1043- 728/1132)  38.7 E30-E19 S17-S08
     10  993  717 (S25W06 L=238)  1064    359        1 ( 818/ 689-1178/ 757)  31.0 E04-W16 S30-S18
     11  469  711 (S39E30 L=201)  1562    354        5 ( 331/ 690- 612/ 739)  15.6 E43-E18 S42-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