Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1089 1619 (N42W20 L=253)  1441    294        3 ( 950/1582-1217/1647)  -7.8 W08-W30 N38-N45
      2  487 1497 (N46E34 L=198)  2053    369        3 ( 342/1431- 658/1571)   6.7 E51-E17 N43-N48
      3 1668 1454 (N17W55 L=288)   773    142        0 (1608/1410-1727/1491) -49.4 W51-W59 N13-N20
      4  249  946 (N14E56 L=177)   665     82        1 ( 230/ 915- 275/ 990) -59.0 E58-E53 N11-N16
      5  537  949 (N10E32 L=201)   317     82        0 ( 509/ 910- 560/ 989) -78.7 E33-E30 N08-N13
      6 1194  727 (S14W05 L=237)   284     82        0 (1178/ 684-1220/ 767)  59.0 W03-W06 S16-S11
      7  371  666 (S05E51 L=182)   493    129        1 ( 355/ 602- 391/ 724)  76.0 E52-E48 S08-N00
      8 1256  603 (S23W06 L=239)   581    147        1 (1203/ 565-1326/ 648) -49.4 W03-W10 S26-S19
      9  613  497 (S18E38 L=195)   575    123        1 ( 575/ 437- 633/ 555)  45.0 E42-E34 S22-S14
     10 1490  413 (S44W22 L=255)  3045    515        5 (1278/ 313-1699/ 530)   7.1 W02-W44 S47-S39
     11  880  338 (S34E26 L=207)   267     74        0 ( 858/ 312- 899/ 377)  33.7 E28-E22 S35-S31
  

    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