Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1660 1649 (N30W72 L= 54)   836     98        1 (1607/1630-1701/1667) -18.4 W65-W77 N28-N32
      2 1374 1538 (N31W35 L= 18)   417     54        1 (1342/1514-1408/1558) -11.3 W33-W38 N29-N32
      3  644 1336 (N31E21 L=321)  1943    473        5 ( 443/1247- 841/1439) -29.4 E39-E07 N22-N40
      4 1108 1129 (N11W07 L=349)   189     54        1 (1086/1103-1129/1160)  33.7 W05-W08 N10-N13
      5  438 1035 (N15E38 L=304)   388    169        0 ( 400/ 984- 509/1070)  18.4 E41-E32 N12-N16
      6  861 1012 (N08E10 L=332)   243     83        1 ( 844/ 969- 882/1050)  68.2 E11-E08 N05-N10
      7  272  942 (N11E54 L=289)  1583    291        5 ( 191/ 830- 346/1048) -45.0 E62-E48 N03-N17
      8 1519  942 (S07W29 L= 12)   278     95        1 (1505/ 908-1547/ 992)  78.7 W28-W30 S08-S03
      9 1161  780 (S10W04 L=347)   778    165        2 (1133/ 697-1187/ 860)  66.0 W01-W06 S14-S05
     10 1384  729 (S17W18 L=360)   450    185        0 (1347/ 696-1447/ 768) -26.6 W15-W21 S19-S14
     11  645  599 (S14E31 L=311)   813    308        0 ( 591/ 540- 717/ 657) -82.9 E34-E26 S17-S09
     12  868  396 (S31E23 L=320)  1178    252        3 ( 784/ 333- 949/ 504)  29.4 E31-E14 S35-S24
  

    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